Ferroplasmas: Magnetic Dust Dynamics in a Conducting Fluid 
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We consider a dusty plasma where the dust particles have a magnetic dipole moment. A Hall- 
MHD type of model, generalized to account for the intrinsic magnetization, is derived. The model 
is shown to be energy conserving, and the energy density and flux is derived. The general dispersion 
relation is then derived, and we show that kinetic Alfven waves exhibit an instability for a low 
temperature and high density plasma. We discuss the implication of our results. 

PACS numbers: 52.25.Dg, 52.25.Xz, 52.27.Lw 



I. INTRODUCTION 



Plasmas with impurity particles, so called dusty or 
complex plasmas, have a manifold of applications, rang- 
ing from applied to basic research problems. Due to this 
wide range of applicability, there has been a steadily in- 
creasing interest in the physics of dusty plasmas (see Ref. 
[l[ and references therein). The above mentioned spec- 
trum of applications covers, e.g., planetary rings, light- 
ning discharges in smoke contaminated air, fusion plas- 
mas, low-temperature laboratory plasmas, and process- 
ing plasmas in the semi-conductor industry. To be spe- 
cific, by a dusty plasma one usually means a three compo- 
nent plasma consisting of electrons, ions, and dust which 
is considered to be significantly heavier than the ions [2| • 
The charge of the dust is assumed to range from a few 
electron charges to thousands. For astronomical applica- 
tions it is often important to also include the dynamics of 
neutral particles [3| . Dusty plasmas contains novel phys- 
ical phenomena, such as dust acoustic waves m, dust 
ion- acoustic waves S, dusty plasma crystals [fj[7j], and 
dust lattice waves [8|, all of which have been experimen- 
tally verified, see Refs. [l|, [9[ and references therein. In 
fact, due to the (in general) relatively low phase veloc- 
ity of dusty plasma waves, these plasmas are useful for 
probing basic properties of plasma excitations. 

Lately there has also been an increasing interest in 
quantum plasma physics [hJ El, E2, EH 0, EH EH , m 
particular the nonlinear aspects of such systems. Such 
plasmas are in general typical for condensed matter en- 
vironments, where the density of the electron gas is high 
[IH E3]j gi vm g a considerable influence from the wave 
function structure of the electrons. Moreover, the spin 
properties of plasmas has been investigated rece ntly by 
means of quantum hydrodynamical models llfj . Il8l 1 1 Oil 
and also by the use of spin kinetic models [20, 123, [221 • 
Even in a regime considered as classical, the effects of 
spin may give a nontrivial influence on the dynamics of 
an electron plasma [23|. Furthermore, spin and the in- 
trinsic magnetic moment of the constituent particles are 
an essential part of magnetic fluids or ferrofluids [241 • A 
fcrrofluid is a mixture of nanosized magnetic particles 
suspended in a liquid. In this article we will consider a 
model where magnetic dust particles are suspended in an 
electron-ion plasma, which can be said to be the plasma 



analogue of a ferrofluid. Such systems has recently been 
investigated both theoretically, considering single par- 
ticle dynamics [H, [2I, [U HH [29| and experimentally 
[3fl l3lf . Here we will consider sufficiently low frequency 
phenomena so that we may use a hydrodynamical model 
to describe the dynamics of the dust particles. 

In Section |TI] we present the governing equation and 
show that the model satisfies an energy conservation law 
in which magnetization transport is included. We then 
linearize the equations in Section IIIII to obtain the gen- 
eral dispersion relation. For the case of a static external 
magnetic field we consider modes propagating perpendic- 
ular to the magnetic field as well as kinetic Alfven waves. 
We show that these modes exhibit instabilities. Finally, 
in Section [TVl we summarize and draw our conclusions. 



II. A MAGNETIZED DUST MODEL 

We here consider a three component plasma consisting 
of electrons, positive ions and negatively charged dust 
particles (denoted by subscripts e, i and d respectively). 
The dust is assumed to be magnetized and have a charge 
—Zde (where e is the elementary charge). The magneti- 
zation of the dust can be assumed to be due to quantum 
mechanical spin, or from a macroscopic magnetization of 
the dust grains themselves. The size of the dust parti- 
cles in, e.g., astronomical environment ranges from a few 
nanometers to about 100 microns and the weight from 
about 10~ 15 g to 10~ 5 g [HQ. In the laboratory for ex- 
ample Ref. [3l| uses grains with size 4.5/^m and mass 
~ 10~ 14 g. Hence the dust is much heavier than the ions, 
m d > mi. 

In the framework of the multi-fluid theory the dynam- 
ics is governed by the continuity and momentum conser- 
vation equations. The continuity equation is given by 



% + V • (n,v a ) = 0, 
at 



(1) 



where the subscript s denotes the different species (s — 
e,i,d), n s is the number density and v s is the fluid veloc- 
ity of particle type s. The momentum conservation for 
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the dust reads 



m d n d I ^ + v d • Vj v d = -Z d en d (E + v d x B) 

-k B V{T d n d ) + M a VB a ,(2) 

where n d is the number density of the dust, k B is Boltz- 
mann's constant, T d is the dust temperature and M is 
the magnetization. We use Einstein's summation con- 
vention for greek indices. The last term in the equation 
above is usually neglected when considering plasmas, but 
for microsized dust grains in a plasma the mutual mag- 
netic dipole interaction can be of importance for labora- 
tory conditions [3(| ■ We will consider perturbations slow 
compared to the plasma frequencies of the ions and the 
electrons. Neglecting the momentum of these particles 
we obtain 

= -e7i e (E + v e x B) -k B V{T e n e ) (3) 
= Z t en t (E + v, x B) - k B V(Tim), (4) 

where T e u-> is the temperature of the electrons (ions) and 
Zje is the charge of the ions. Note, we neglect the spin of 
the electron and the protons. The magnetization of the 
dust is assumed to be orders of magnitude larger. The 
equations above are coupled to Maxwell's equations 



that (T e +Ti/Zi)n e <C (T d +Z d Ti/ 1 Zi)n d in order to obtain 
a closed set of equations. This condition can of course 
be relaxed, introducing new thermal effects. However, 
the principle dynamics of the complex plasma will not 
be significantly affected. Using some vector identities we 
can finally write the continuity and momentum equations 
for the dust 



^ + V(n d v d ) = 0, 



(11a) 



and 



m d n d 



-V 



d 



v d • V v d = B • V 



M 

Mo 

MB-:- k L; | L; + ^^jn d 



(lib) 



respectively. Using Eqs. (HJ), and (JSJ) we can derive 
the time evolution equation for the magnetic field 



<9B 

~dt 



V x (v d x B) + V x 



[V x (B - noM)] x B 

(j, Z d en d 



(11c) 

The magnetization is taken to be proportional to the den- 
sity of the dust particles and in the direction of the mag- 
netic field 



V-E = - — (n e - ZiUi + Z d n d ), (5) 
VxE = --, (6) 



M = fi d (B,T d )n d B. 



(lid) 



and 



V x B = (LtoV x M - e[i (n e v e - ^r^Vi + Z d n d v d ). (7) 

A closed system of equations for the dust can be derived. 
To do this we start by adding the momentum Eqs. ([3|) 
and (j4]) and assume that quasi- neutrality holds Z^rii rs 
n e + Z d n d to obtain 

-eZ d n d ~E = -e(n e v & — ZiUiVi) x B - fc B V(Tjni + T e n e ). 

(8) 

Solving this for the electric field and inserting it into the 
momentum equation for the dust @ yields 

( d w\ 
m d n d \—+v d -\7jv d — 

-e(n e v e - Z i n i v l + Z d n d v d ) x B 
~k B V{T d n d + T t n t + T e n e ) + M a VB a . (9) 



Using Eq. ([7]) we obtain 
d 



m d n d 



-k R y 



Of 



v d • V )v d = V x 



B 

Mo 



T, 



Ti 
Z, 



T d +—— I n d 



M ) x B (10) 
M a VB a 



where the quasi-neutrality condition has been used to 
rewrite the thermal pressure terms and we will assume 



This model includes the case where the particles have an 
intrinsic magnetic moment that will be aligned with an 
applied magnetic field, which is what we have in mind. 
However, it also includes the case where the dust parti- 
cles have no net magnetic moment. The occurrence of an 
external magnetic field will induce a magnetization in the 
dust particles which yields a macroscopic magnetization 
of the fluid. A magnetization of the dust could also arise 
from spinning of the dust particles (with charges attached 
to the surface) which yields a diamagnetic response to an 
applied field. However, the magnetic dust-dust interac- 
tion due to this can often be neglected. See Ref. [28| for 
a more detailed discussion. Our model can be compared 
with a ferrofluid which is a colloidal suspension of mag- 
netic particles in a liquid. In an ionic ferrofluid the mag- 
netic particles are kept apart by repulsive electrostatic 
forces, see e.g. Refs. [3, [3l|, |H, 0, [37j • Due to the high 
density and the correspondingly high collision frequency 
we note that free currents can typically be neglected in 
ionic ferrofluids. By contrast, the simultaneous existence 
of free currents and magnetic dipole moment have been 
shown to be significant in dusty plasmas [311 ]. A more 
general assumption than Eq. (jlldp would be to assume 
that the magnetization is also dependent on the electron 
and ion temperatures since collisions with these particles 
may change the magnetization. 

With the exception for the occurrence of the magneti- 
zation due to the magnetic moment of the dust, M, these 
equations are the same as Hall-MHD theory [32|, [HI]. 
Moreover, it should be noted that the structure of the 
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system of Eqs. (lla)-(llc) is the same as one get from 
the magnetized ideal MHD model of Refs. jiol |34| . 
if that system is extended to include the Hall current. 
In that case naturally the dust density and velocity in- 
stead will refer to the ion density and velocity. Although 
Eqs. (lla)-(llc) can describe an electron-ion plasma, it 
should be noted that the physics for that case is different 
in several respects: Firstly, in the case of an electron- 
ion plasma, it is the lighter species, the electrons, that 
contribute to the magnetization, due to their magnetic 
moments being larger than those of the ions. Secondly 
we stress that the validity conditions of the model have 
no simple correspondence between the dust dominated 
plasma and the electron-ion plasma case. In what fol- 
lows, we will mainly be concerned with the dusty plasma 
applications. 

Next we want to find an energy conservation law for 
the Eqs. (TTTj) . In order to do this we must specify an 
equation of state for the system. For simplicity we choose 
the simple model 



P 



ndo 



(12) 



is the energy per volume and 



p = "*^3 V + jg_ v _ (B . M )v 



2 

v x B 



7-1 
1 

Zden d 



(V x H) x B 



x H (15) 



is the flow of energy out of the region. In Eq. (|14p the 
first term is the kinetic energy per volume, the second 
term is the energy density from the pressure, the third 
term is the energy stored in the magnetic field and the 
last term is the energy in each volume element due to the 
magnetic moment of the dust particles. Similarly, in Eq. 
(|15[) . the first three terms are the flow of kinetic, pressure 
and magnetic energy density that follows the flow of each 
volume element. The last term is the Poynting vector 
which is modified by the inclusion of the Hall-term and 
the magnetization. 



for the total pressure P = kg(Td + Z^Ti/Zi), where Po 
and rido are the equilibrium pressure and density respec- 
tively. With this equation of state the energy conserva- 
tion becomes 



III. THE DISPERSION RELATION 



where 



W 



dW 
~dt 



m d n d v^ 



V-P = 0, 



P 



B*_ 

2fiQ 



B M 



(13) ii near j ze Eq S (jTTJ) and Fourier decompose. Fur- 
thermore, the coordinate system is defined so that Bo = 
BqZ and k = fc^x + k z z. For simplicity we choose an 
isothermal pressure model. This gives the dispersion re- 

(14) lation 



kUVl-V d %) 



• da 

i 

&cd 



-i—(k 2 v^ A -k 2 X B ) 

Ucd 



k 2 V 2 
K - v dA 



k 2 V 2 
h z v dA 



-k x k z V^ a 



i k x k z VdM 

&cd 



k x k z V dc 



U c d 



k x k z VdA 



= 0, 



(16) 



-k 2 V 



z v da 



where V 2 A = V 2 A - V 2 M , V 2 a = V 2 a - V 2 M and we have 
defined 



o kn (' ZdT« 
m d \ Zi 



V? 



V dA = 



V 2 
v dM 



• 113 



Ho m dn do 
fJ-doBo 

m d 
dfJ-do Bq 
dB m d 



(17a) 
(17b) 
(17c) 
(17d) 



The velocity Vd a is a generalized thermal speed for the 
dust, V d A is the dust Alfven speed and VdM and VdB are 
related to the magnetization of the dust. The frequency, 
uj c d = ZdeBo/md is the cyclotron frequency of the dust. 
The dispersion relation, Eq. (|16p . is in general a third 
degree polynomial in ui 2 . Specifically, for oj <C u> c d the 
three roots to the dispersion relation are the fast and slow 
magnetosonic modes, and the shear Alfven wave [38|. 

To see the implications of the derived dispersion rela- 
tion a couple of special cases are now considered. For 
a wave propagating perpendicular to the magnetic field 
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(k = fc^x) the dispersion relation is obtained from Eq. 
(HHJ) and reads 



= kl Vf A + V r . 



da 



v. 



dB 



(18) 



To get a qualitative description of the instability con- 
dition for Eq. (fT8"|) we assume that the dust and the 
ions are in thermal equilibrium so that we may assume 
(T d + Z d T t /Zi) = NT d where N > 1 is a constant. Typ- 
ically we have Z d jZi ranging from unity to a few thou- 
sands Furthermore we assume that the spins are 
thermally distributed (see e.g. Ref. [38j) the magnetiza- 
tion per density can be written 



= Hd tanh 



( fi d B\ 
\k B T d ) 



(19) 



In the case that the dust particles have a high total 
spin number the tanh-function in Eq. (|19p correspond- 
ing to Fermi-Dirac statistics should be replaced by the 
Langevin-function, corresponding to Maxwell-Boltzmann 
statistics. The difference between these functions are rel- 
atively small, however, and thus we will use Eq. (fT9|) for 
the remainder of this paper. We note that the magnetic 
moment fi d can be several orders of magnitude larger 
than the Bohr magneton. The condition for instability, 
Eq. |[TB]). then becomes 



B 



N 



k B Td u d Q n d B a 



1 



fid 



<0, 



Hmdol+d HdB fi d k B T d 

J (20) 

where //<jo = Mrf(-^o) m accordance with Eq. (p~5|) . Note 
that this equation implies that more of the magnetic 
dipoles are re-oriented towards the lower energy state 
when the magnetic field strength is increased. For this 
to apply, the relaxation time to reach the lower energy 
state must be shorter than the wave period time. In case 
the opposite ordering holds, the fraction of particles in 
the different energy states remain constant during a wave 
period, and consequently the term proportional to V^ B 
in Eqs. ([To]) and (fTg)) should be dropped, which corre- 
spond to neglecting the fourth and last term in Eq. (|2U)l . 
The difference in the dispersion relation, depending on 
whether the magnetic dipoles have time to change dur- 
ing a wave period or not, turns out to be relatively small, 
however. For defmitcness, we will stick to the case where 
the relaxation time of the magnetic dipoles is sufficiently 
fast for Eq. (f2"0"]) to apply for the rest of this work. 

Next we consider the kinetic Alfven type of waves. 
In this case the ordering k z <C k x , V da <C V d A and 
^ ~ k z V d A applies. The dispersion relation can then 
be approximated by 



(21) 



,2 kl 



1+- 



dA~ 



v^ B )vl + vi M (vl 



'dMVdA 



-V? - 

da 



'dB 



(VdA + V da 



For V? M , V dB — > this reduces to the well known kinetic 



Alfven waves [39| . Analyzing Eq. (|2"Tj) it is seen that the 
wave-mode can be unstable provided the numerator of 
the second term of the right hand side is negative, i.e. 
we obtain the instability condition 

{v!a - visWl + v^ M (vi A + vl - vi B ) < o (22) 

Assuming once more that the dust and the ions are in 
thermal equilibrium we get 



B, 



o 



fi d B , (N-l)fi d B 



N k B T d 



(23) 

The conditions for instabilities Eqs. (|2"D|) and (|2"5)) have 
been plotted in Fig. [TJ Note that in order to have insta- 
bilities we need to have sufficiently high densities and/or 
sufficiently low temperatures. 

We here give the following simplified picture of why the 
instability of this type can occur. The volume elements of 
the plasma are electrically neutral since the electron and 
ion background will screen any excess electrical charge. 
Further, the magnetization of a volume element is in the 
direction of the magnetic field and the different volume 
elements will attract each other like small magnets. Con- 
sider now the magnetic flux through a surface with nor- 
mal parallel to B. If the oscillations have k x = then 
the magnetic flux through the surface will not change, 
and hence there will be no build up of magnetization. If, 
on the other hand, the oscillations occur perpendicular 
to the field there can be a local build up of the magnetic 
field. This can for sufficiently low temperature and high 
density cause the plasma to collapse similarly to the case 
of the Jeans instability [38l . 140 ] . 



IV. SUMMARY AND CONCLUSION 

In the present paper we have put forward a Hall MHD 
type of model with an intrinsic magnetization. We have 
shown that the magnetized Hall MHD model is energy 
conserving, and presented the expressions for energy den- 
sity and the energy flux. A set of equations of this 
type could describe different types of systems: An or- 
dinary electron-ion plasma, in which case the magnetiza- 
tion would be due to the electron spin, or - as emphasized 
here - a three component plasma containing electrons, 
ions and heavy dust particles. 

As the next step, we have investigated the linear modes 
and the stability properties of the homogeneous system. 
The general dispersion relation has been derived, describ- 
ing the fast and slow magnetosonic modes and the shear 
Alfven waves, as modified by the Hall current and the 
magnetization of the system. Due to the magnetization, 
the homogeneous system may be unstable, as predicted 
already from an magnetized ideal MHD type of model 
[38| . The main new finding from the stability analysis in 
this paper, is that inclusion of the effects due to the Hall 
current extends the unstable region of parameter space, 
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(a) 



(b) 



Figure 1: The instability regions for the (a) the perpendicular mode, Eq. (|20p and the (b) kinetic Alfven mode, Eq. (|23|) plotted 
in terms of the dimensionless parameters X = ftdBo / (fesTd) and Y = Bo/(fJ>ondop>d) with N = 10. 



as described by Fig 1 . The wave that first becomes unsta- 
ble turns out to be the magnetized version of the kinetic 
Alfven wave. 

For the effects of magnetization to be significant, we 
need relatively high density plasmas, and/or low temper- 
atures. For the case of electron-ion plasmas, these can be 
found in astrophysical environments, such as the interior 
of white dwarf stars and pulsars. 

We can make some numerical estimations for the pa- 
rameters X = p, d B /(k B T d ) and Y = B / (p n do p. d ) . 
From Ref. [3C| we find that a the magnetic moment per 
particle can be of the order 10 _12 m. _2 A _1 and magnetic 
induction of the order 0.1T. Furthermore, we assume 
that the dust temperature is low T d ~ IK. The de nsity 
of particles is taken to be n d o ~ 10 12 m -3 as in Ref. [261 ]. 
This gives us X ss 10 10 and Y w 10 6 . Comparing this 
with Figs. [T] we see that the instabilities considered here 
are not possible to detect in current experiments. As 
shown e.g. by Refs. [3(| HH, however, magnetic dipole 



effects can still be of significance in dusty plasmas. 

The model developed here should be considered only 
as a first step since it does not account for potentially im- 
portant effects of more elaborate models. These include 
the two-fluid model, where spin up and spin down popu- 
lations are considered as different species [23| , the kinetic 
description 12211 and models including nearest neighbor 
interactions [29l l3fj| , which is important in the strongly 
coupled regime. 
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